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LONG-TERM  GOALS 

The  long  term  goals  of  this  projeet  are  to  develop  stoehastie  models  for  broadband  mode  signals 
propagating  in  fluetuating  oeean  environments  and  to  design  robust  signal  proeessing  techniques  for 
these  ocean  environments. 

OBJECTIVES 

To  develop  a  greater  understanding  of  broadband  mode  signals  in  fluctuating  ocean  environments,  this 
project  focused  on  several  closely -related  research  objectives.  The  first  objective  was  to  characterize 
mode  fluctuations  at  megameter  ranges  using  experimental  measurements  and  simulations.  Mode 
arrivals  at  long  range  are  very  complicated  due  to  intemal-wave-induced  coupling.  Erom  a  theoretical 
standpoint,  the  effects  of  internal  waves  on  long-range  sound  propagation  are  not  thoroughly 
understood,  and  much  of  the  previous  work  in  this  area  focused  on  the  ray  arrivals  because  they  are 
amenable  to  analysis  via  the  geometrical  optics  approximation.  Characterizing  the  fluctuations  is  the 
first  step  towards  the  long-term  goal  of  developing  a  sufficiently  general  model  of  mode  propagation 
that  explains  recent  experimental  observations  and  clarifies  the  relationship  between  mode  and  ray 
representations  of  the  field.  The  second  objective  of  this  project  was  to  develop  a  framework  for  mode 
processing  that  mitigates  the  effects  of  environmental  mismatch,  sensor  failures,  and  interference.  Since 
mismatch  is  a  problem  that  plagues  many  sonar  applications,  this  work  has  implications  beyond  the 
current  project. 

APPROACH 

To  address  the  stated  objectives,  this  project  combined  experimental  data  analysis  with  simulation  and 
theory.  Over  the  last  10  years,  a  series  of  ONR- sponsored  experiments  generated  rich  data  sets  for 
studying  normal  mode  signals  at  megameter  ranges.  The  Acoustic  Thermometry  of  Ocean  Climate 
(ATOC),  Alternate  Source  Test  (AST),  and  North  Pacific  Acoustic  Eaboratory  (NPAE)  experiments 
were  the  primary  focus  of  this  project.  The  ATOC,  AST,  and  NPAE  data  sets  include  vertical  array 
receptions  at  center  frequencies  of  28  Hz,  75  Hz,  and  84  Hz  at  ranges  of  3.5  Mm  and  5.1  Mm.  As  a  part 
of  this  project  the  PI  and  her  students  analyzed  these  data  sets  using  signal  processing  techniques 
developed  for  processing  the  low  order  modes.  The  PI  also  contributed  to  the  planning  and 
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implementation  of  the  2004-2005  SPICE04/LOAPEX  experiment  that  finished  in  June  2005. 


Between  2002  and  2005,  several  George  Mason  University  eleetrieal  engineering  students  worked  on 
this  projeet.  Mr.  Tarun  Chandrayadula  investigated  the  problem  of  designing  adaptive  equalizers  for 
broadband  mode  signals.  Mr.  Aravinda  Sringarapuram  worked  on  mode  estimation  methods  that  ean 
adapt  to  the  ehanging  interferenee  eharaeteristies  of  the  ehannel.  ET  Erie  St.  Pierre  attended  GMU 
sponsored  by  the  U.S.  Coast  Guard  postgraduate  training  program.  In  2004-2005,  he  developed  a 
method  for  estimating  the  sound  speed  profile  required  for  mode  proeessing  from  sparse  measurements 
of  the  temperature  along  a  vertieal  array. 

WORK  COMPLETED 

The  duration  of  this  projeet  was  3.5  years,  from  April  2002  to  Deeember  2005.  This  report  summarizes 
the  work  eompleted  and  the  publieations  produeed  during  the  projeet. 

ATOC/AST  Mode  Analysis:  An  extensive  statistieal  analysis  of  the  mode  arrivals  in  the  ATOC  and  AST 
experiments  was  eompleted  in  2003-2004.  Two  artieles  by  Wage  et  al.  in  the  Journal  of  the  Acoustical 
Society  of  America  deseribe  the  mode  statisties,  sueh  as  eoherenee  and  time  spread,  eomputed  from 
these  experiments  [1,2].  Kathleen  Wage  summarized  these  results  in  an  invited  talk  at  the  November 
2004  meeting  of  the  Aeoustieal  Soeiety  of  Ameriea  [3].  She  also  presented  interim  reports  on  this  work 
at  three  North  Paeifie  Aeoustie  Eaboratory  Workshops  between  2002  and  2004  [4,  5,  6]. 

Approximate  Mode  Filtering:  Mode  filters  are  typieally  designed  using  modeshapes  derived  from 
numerieal  solution  of  the  Helmholtz  equation  for  a  measured  sound  speed  profile  (SSP).  In  experiments 
sueh  as  NPAE,  the  temperature  (henee  the  sound  speed)  is  only  measured  aeross  the  span  of  the  array. 
To  faeilitate  mode  proeessing,  it  is  helpful  to  have  an  approximate  solution  for  the  modes  that  only 
requires  the  environmental  information  along  the  array.  Kathleen  Wage  presented  a  paper  at  the  2004 
Asilomar  Conferenee  on  Signals,  Systems,  and  Computers  that  explored  the  use  of  a  uniform  WKB-like 
approximation  to  design  mode  filters  [7].  She  presented  the  initial  results  of  this  work  at  the  Underwater 
Aeoustie  Signal  Proeessing  Workshop  in  2003  [8]. 

Sound  Speed  Estimation  from  Sparse  Environmental  Measurements:  Mode  filtering  requires  aeeurate 
knowledge  of  the  sound  speed  profile  at  the  array.  Eor  long-duration  experiments,  the  environment  is 
typieally  measured  using  a  relatively  small  set  of  temperature  sensors  loeated  along  the  array.  In  the 
2004-2005  a  method  was  developed  for  extrapolating  these  sparse  temperature  measurements  to  obtain 
a  sound  speed  profile  for  use  in  mode  proeessing.  The  approaeh  was  applied  to  temperature 
measurements  aequired  during  the  NPAE  experiment.  Erie  St.  Pierre  presented  a  paper  deseribing  these 
results  at  the  September  2005  lEEE/MTS  Oeeans  Conferenee  [9]. 

Mode  Equalization:  At  megameter  ranges,  the  low  mode  arrivals  are  quite  eomplieated.  As  a  part  of 
this  projeet,  the  problem  of  equalizing  these  signals  has  been  studied.  Tarun  Chandrayadula  presented  a 
paper  on  the  design  of  mode  equalizers  at  the  September  2005  lEEE/MTS  Oeeans  Conferenee  [10].  He 
also  presented  a  talk  on  the  multi-mode  deeision  feedbaek  equalizer  for  underwater  aeoustie 
eommunieations  at  the  2005  Underwater  Aeoustie  Signal  Proeessing  Workshop  [11]. 

SPICE04/LOAPEX  Experiment:  PI  Kathleen  Wage  eontributed  to  planning  diseussions  for  the  2004 


experiment  at  several  NPAL  workshops.  In  particular  she  analyzed  the  mode-resolving  capabilities  of 
the  proposed  array  geometries.  These  results  were  presented  at  NPAL  workshops  in  2002  [12]  and 
2003  [13],  and  summarized  in  a  memo  written  in  April  2004  [14].  Kathleen  was  also  a  member  of  the 
science  team  for  the  SPICE04  deployment  and  recovery  cruises,  which  took  place  in  June  2004  and 
June  2005,  respectively.  The  chief  scientist  for  these  cruises  was  Peter  Worcester  (Scripps  Institution  of 
Oceanography).  Two  vertical  line  arrays  (VLAs)  and  two  source  moorings  were  deployed  for  the 
SPICE04/LOAPEX  experiment. 

RESULTS 

The  following  paragraphs  summarize  key  results  obtained  from  2002  to  2005.  Parts  of  this  summary 
have  appeared  in  previous  annual  reports. 

ATOC/AST  Mode  Analysis:  A  key  contribution  of  this  project  was  the  statistical  analysis  of  the  mode 
arrivals  in  the  ATOC  and  AST  data  sets.  This  analysis  is  thoroughly  described  in  two  journal  articles  by 
Wage  et  al.  [1,2].  Eigures  1  and  2  are  representative  of  the  mode  statistics  described  in  these  two 
papers.  Eigure  1  shows  the  time  spread  measured  at  the  Hawaii  and  Kiritimati  VEA’s,  which  were 
located  at  ranges  of  3.5  Mm  and  5.1  Mm,  respectively.  The  figure  illustrates  increase  in  time  spread  as  a 
function  of  mode  number.  The  28  Hz  mode  1  measured  at  Hawaii  show  the  least  amount  of  spread.  At 
28  Hz  the  internal  wave  scattering  is  significantly  reduced,  and  mode  1  is  dominated  by  a  single  arrival. 
At  higher  frequencies  and  longer  ranges,  internal  wave  scattering  produces  a  more  complicated  arrival 
pattern.  Note  that  the  average  time  spread  for  the  first  10  modes  is  approximately  1.5  seconds,  which 
corresponds  to  a  coherence  bandwidth  of  0.67  Hz  for  the  channel.  Eigure  2  compares  the  temporal 
coherence  of  mode  1  measured  using  receptions  from  the  75  Hz  ATOC  source  and  the  28  Hz  AST 
source.  The  coherence  time  is  significantly  longer  for  the  28  Hz  receptions  than  for  the  75  Hz 
receptions.  Both  show  reasonable  agreement  with  simulations  of  propagation  through  internal  waves. 
See  the  paper  by  Wage  et  al.  for  more  discussion  of  temporal  coherence  and  how  to  measure  it  [2]. 

Approximate  Mode  Filtering:  Spatial  filtering  for  modes  requires  knowledge  of  the  modeshapes,  which 
requires  accurate  knowledge  of  the  environment.  Unfortunately,  environmental  information  is  often 
incomplete.  Eor  example  in  both  the  1998  North  Pacific  Acoustic  Laboratory  (NPAL)  experiment  and 
the  SPICE04  experiment,  the  environment  was  sampled  along  the  VLA  span,  but  not  above  and  below 
the  array.  To  facilitate  mode  processing,  an  approximation  to  the  modeshapes  is  needed  that  only 
requires  the  environmental  information  along  the  array.  Ideally,  the  approximation  should  provide  an 
analytical  form  for  the  modeshapes,  rather  than  just  a  numerical  solution.  An  analytical  form  is 
desirable  because  it  provides  more  intuition  about  the  filtering  problem  and  may  permit  simpler  design 
of  robust  mode  filters.  WKB  theory  provides  the  most  common  approximate  solution  to  the  mode 
problem.  The  main  problem  with  using  the  standard  WKB  solution  to  design  mode  filters  is  that  it  is 
singular  at  the  turning  points.  Miller  and  Good  developed  a  uniform  WKB-like  solution  to  the 
Schrodinger  equation  that  is  continuous  at  the  turning  points  [15].  The  paper  presented  by  Wage  at  the 
Asilomar  Conference  in  2004  describes  the  application  of  Miller  and  Good’s  method  to  the  underwater 
acoustic  mode  problem  [7].  The  major  finding  of  this  paper  is  that  the  approximate  mode  filter  may 
perform  better  (in  terms  of  crosstalk)  than  the  exact  mode  filter  when  the  approximate  filter  can  be 
designed  using  measured  sound  speed  across  the  VLA  and  the  exact  mode  filter  can  only  be  designed 
using  an  archival  sound  speed  measurement. 


Estimated  time  spread 


Figure  1:  Time  spread  for  the  first  10  modes  measured  at  the  Hawaii  and  Kiritimati  VLAs,  denoted 
HVLA  and  KVLA,  respectively.  The  75  Hz  results  are  for  the  bottom-mounted  Pioneer  Seamount 
source  deployed  for  the  ATOC  experiment.  The  28  Hz  results  are  for  the  ship-suspended  source 
deployed  during  the  AST  experiment.  The  average  time  spread  is  1.5  seconds,  and  the  spread 
increases  with  mode  number.  For  a  complete  discussion  of  the  mode  time  spread  measured  during 
the  ATOC  and  AST  experiments,  see  the  paper  by  Wage  et  al.  [2]. 


MSC  for  mode  1 


Figure  2:  Comparison  of  magnitude-squared  coherence  for  mode  1  at  75  Hz  and  28  Hz  measured 
during  the  ATOC  and  AST  experiments  with  simulation  results.  The  28  Hz  data  show  significantly 
longer  coherence  times  than  the  75  Hz  data,  and  there  is  good  agreement  with  simulation  of 
propagation  through  internal  waves.  The  paper  by  Wage  et  al.  describes  these  results  in  more 

detail  [2]. 


Sound  Speed  Estimation  for  Mode  Processing:  As  noted  above,  mode  filtering  of  underwater  aeoustie 
reeeptions  depends  on  aeeurate  knowledge  of  the  mode  shapes  and  wavenumbers.  From  an 
experimental  standpoint,  it  is  often  impraetieal  to  measure  the  environmental  parameters  that  determine 
the  SSP  on  a  dense  grid  in  both  time  and  depth.  For  example,  in  long-duration  aeoustie  experiments 
sueh  as  NPAL,  high  resolution  measurements  of  temperature  and  salinity  are  only  possible  during  the 
deployment  and  reeovery  eruises.  In  the  intervening  months  the  environment  is  sampled  by  a  relatively 
small  number  of  sensors  attaehed  to  the  aeoustie  mooring.  One  approaeh  to  defining  the  SSP  for  mode 
filter  design  is  to  interpolate  and  extrapolate  these  environmental  measurements.  Other  approaehes 
would  be  to  use  arehival  environmental  data  or  the  output  of  an  oeean  model.  The  main  problem  with 
arehival  data  and  oeean  model  data  is  that  it  ean  be  severely  biased.  In  2004-2005,  we  developed  a 
method  for  interpolating  and  extrapolating  the  NPAL  temperature  measurements  to  obtain  sound  speed 
profiles.  The  teehnique,  deseribed  in  [9],  assumes  that  the  ehanges  in  the  SSP  assoeiated  with 
mesoseale  and  seasonal  effeets  ean  be  written  in  terms  of  the  quasi-geostrophie  dynamie  modes.  The 
dynamie  modes  are  an  orthonormal  basis  derived  from  historieal  oeeanographie  data  for  the  VLA 
loeation.  By  fitting  the  observed  temperature  fluetuations  to  the  dynamie  modes,  it  is  possible  to  obtain 
an  estimate  of  the  SSP  required  for  mode  proeessing.  The  Spaee-Time  Kalman  Filter  (STKF)  provides 
a  eonvenient  framework  for  solving  the  estimation  problem.  The  paper  by  St.  Pierre  and  Wage  [9] 
deseribes  the  STKF  implementation  and  examines  its  performanee  for  the  NPAL-98  data  set.  It 
analyzes  the  sound  speed  errors  assoeiated  with  the  method  and  eonsiders  how  those  errors  propagate 
through  the  mode  eomputations. 

Mode  Equalization:  Tomographie  experiments  use  broadband  low  frequeney  signals  sueh  as  maximal 
length  sequenees  (M-sequenees)  to  study  the  oeean.  One  problem  in  deteeting  the  mode  arrivals  at 
megameter  range  is  that  multipath  eauses  time  spread  in  the  transmitted  M-sequenee,  leading  to 
interferenee  between  adjaeent  symbols.  This  problem,  ealled  Inter  Symbol  Interferenee  (ISI),  is 
eommon  in  wireless  ehannels.  Equalizers  are  used  in  underwater  aeoustie  eommunieations  to  reduee 
ISI  and  provide  SNR  gain  e.g.,  see  [16].  Working  with  data  from  the  ATOC  Engineering  Test,  Ereitag 
and  Stojanovie  demonstrated  that  deeision  feedbaek  equalizer  (DEE)  ean  be  used  for  very  long  range 
underwater  aeoustie  signals  [17].  During  the  ATOC  Engineering  Test,  M-sequenee  signals  were 
reeorded  aeross  a  20-element  VEA.  Ereitag  and  Stojanovie’s  equalization  seheme  eonsisted  of  a 
multi-ehannel  DEE  that  operated  on  the  multiple  ehannels  of  the  VEA  and  sueeessfully  deeoded  the 
symbols  with  zero  deeoding  error. 

Motivated  by  this  earlier  work,  we  began  exploring  the  design  of  equalizers  for  the  low  mode  signals.  In 
2004-2005,  we  implemented  a  DEE-based  equalizer  for  the  low  modes  [10].  In  tomography  (as  opposed 
to  eommunieations),  the  transmitted  symbols  are  known.  Eor  the  mode  equalizer  the  transmitted 
symbols  are  fed  baek  instead  of  the  symbol  deeisions,  henee  the  mode  equalizer  is  equivalent  to  the 
DEE  operating  in  training  mode.  The  output  SNR  of  the  equalizer  provides  a  useful  metrie  for 
analyzing  the  mode  arrivals.  The  mode  equalizer  works  as  follows.  Eirst,  a  spatial  filter  is  used  to 
estimate  the  desired  mode  time  series  that  the  equalizer  will  operate  on.  Then  the  equalizer  is  allowed  to 
synehronize  to  a  series  of  time  instants  and  the  output  SNR  is  ealeulated  for  eaeh.  When  the  equalizer 
synehronizes  to  a  time  instant  where  it  is  unable  to  suoeessfully  eombine  multipaths,  it  will  deeode 
symbols  with  a  high  mean  square  error,  thus  the  output  SNR  will  be  low.  When  the  equalizer 
synehronizes  to  a  time  instant  that  enables  it  to  eombine  all  the  multipaths  in  the  signals,  the  deeoded 
symbols  will  have  a  low  mean  square  error,  thus  a  high  SNR.  Chandrayadula  presented  a  paper  at  the 
2005  lEEE/MTS  Oeeans  Conferenee  [10]  demonstrating  that  equalizers  ean  be  designed  for  the  low 


mode  signals  in  the  NPAL  receptions.  Figure  3  shows  examples  of  SNR  vs.  synchronization  time 
curves  for  modes  1  and  10  in  one  of  the  NPAL-98  receptions  (for  the  40-element  VLA  located  off 
California).  The  plots  illustrate  how  the  results  vary  when  equalizers  of  different  lengths  are  used. 
Ideally  the  equalizer  length  should  be  matched  to  the  multipath  spread  of  the  modes.  Since  the  SNR 
keeps  increasing  for  equalizers  up  to  2  seconds  long,  it  appears  that  the  multipath  spread  of  the  modes 
in  the  NPAL-98  receptions  is  on  the  order  of  2  seconds  or  greater. 

Long-Range  Communications:  In  addition  to  investigating  how  equalizers  can  be  used  for  tomography, 
this  project  also  developed  a  multi-mode  decision  feedback  equalizer  (MMDFE)  for  communications 
applications  and  applied  it  to  data  from  the  NPAL-98  experiment.  The  MMDFE  is  based  on  the  reduced 
complexity  multi  channel  DFE  by  Stojanovic  et  al.  [18].  Figure  4  shows  the  data  processing  path 
required  for  implementing  the  MMDFE  on  the  NPAL-98  data.  As  shown,  the  VLA  receptions  are  first 
corrected  for  clock  drift  and  mooring  motion.  Then  a  least  squares  mode  beamformer  similar  to  the  one 
designed  by  Wage  [2]  is  applied  to  obtain  the  estimated  time  series  for  the  first  10  modes.  The  mode 
filtering  provides  spatial  array  gain  and  reduces  the  computational  complexity  by  a  factor  of  two.  The 
reduction  in  computational  complexity  is  achieved  in  two  ways.  First,  mode  filtering  reduces  the 
number  of  channels  the  DFE  operates  on.  Second,  mode  filtering  also  facilitates  a  reduction  in  the 
number  of  taps  required  for  the  DFE.  To  see  this,  note  that  the  received  pressure  field  is  a  linear 
combination  of  the  mode  signals.  A  multi-channel  DFE  operating  on  multiple  channels  of  the  VLA  has 
to  contend  with  more  amount  of  multipath  than  a  DFE  operating  on  the  mode  signals.  The  multi-mode 
DFE  can  thus  operate  with  a  fewer  number  of  taps  for  each  channel.  Figure  5  shows  the  scatter  plots  of 
the  equalized  symbols  for  various  configurations  of  the  MMDFE.  The  columns  show  the  results  for 
using  1  mode,  5  modes,  and  10  modes  in  the  equalizer.  The  rows  show  the  results  for  using  different 
amounts  of  signal  averaging  prior  to  implementing  the  MMDFE.  As  Figure  5  shows,  the  DFE  operating 
on  a  single  mode  does  not  have  enough  spatial  diversity  to  separate  the  symbols.  Additional  time 
averaging,  which  should  improve  the  SNR,  does  not  improve  the  performance  of  the  single-mode  DFE. 
In  contrast,  using  20-period  averages  appears  to  substantially  improve  the  performance  of  the  5-mode 
and  10-mode  DFEs.  These  scatter  plots  illustrate  how  the  MMDFE  is  capable  of  exploiting  the  spatial 
diversity.  Figure  6  further  illustrates  the  effect  of  temporal  averaging  on  the  performance  of  the 
MMDFE.  The  plot  shows  the  Bit  Error  Rate  (BER)  versus  the  number  of  M-sequence  periods  averaged. 
The  BER  curves  show  that  the  performance  of  the  10-mode  DFE  does  not  significantly  differ  from  the 
5-mode  DFE.  The  fact  that  the  DFE  benefited  from  temporal  averaging  indicates  that  the  multipath 
structure  for  the  modes  is  coherent  over  multiple  M-sequence  periods. 

IMPACT/APPLICATIONS 

The  objectives  of  this  research  have  both  scientific  and  operational  applications.  The  statistical 
characterization  of  the  acoustic  channel  provides  crucial  information  to  guide  the  design  of  long-range 
systems  for  tomography,  communication,  and  surveillance.  The  signal  processing  techniques  being 
developed  are  applicable  in  both  short-  and  long-range  propagation  scenarios. 

RELATED  PROJECTS 

This  work  is  closely  related  to  the  North  Pacific  Acoustic  Laboratory  project,  directed  by  principal 
investigators  Peter  Worcester  (SIO)  and  James  Mercer  (APL-UW).  A  number  of  other  ONR-sponsored 
researchers  work  on  projects  related  to  NPAL  and  participate  in  the  NPAL  Workshops. 
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(a)  Mode  1  time  series  (4-period  avg). 
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(c)  Mode  10  time  series  (4-period  avg). 
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Figure  3:  Example  of  mode  time  series  measured  in  the  NPAL-98  experiment  and  the  corresponding 
mode  equalizer  outputs.  The  plots  on  the  left  show  the  mode  1  and  10  time  series  estimated  from  a 
4-period  average  of  the  received  signal.  The  plots  on  the  right  show  the  output  SNR 
vs.  synchronization  time  curves  obtained  from  the  mode  equalizer.  See  the  paper  by  Chandrayadula 

and  Wage  for  more  details  about  these  results  [10]. 
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Figure  4:  Block  diagram  of  the  multi-mode  decision  feedback  equalizer  (MMDFE)  for  the  NPAL-98 
experimental  data.  As  shown,  processing  consists  of  mooring/timing  corrections,  followed  by  mode 
beamforming,  and  then  decision  feedback  equalization.  Mode  beaforming  provides  spatial  array  gain 
and  reduces  the  required  dimension  of  the  decision  feedback  equalizer. 
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Figure  5:  Scatter  plots  for  various  configurations  of  the  MMDFE  operating  on  NPAL-98  data.  The 
rows  show  results  for  different  amounts  of  signal  averaging  prior  to  processing.  The  columns  show 
results  for  different  numbers  of  modes  included  in  the  equalizer.  This  result  was  presented  at  the 
2005  Underwater  Acoustic  Signal  Processing  Workshop  [11]. 


BER  performance  of  the  MMDFE 


Figure  6:  Bit  error  rate  as  a  function  of  the  number  of  M-sequence  periods  averaged  prior  to 
processing.  This  result  was  presented  at  the  2005  Underwater  Acoustic  Signal  Processing 

Workshop  [11]. 
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